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A gate electric field has a small but non-negligible effect on the phase of the valley-orbit coupling 
in Si quantum dots. Finite interdot tunneling between valley eigenstates in a double quantum dot is 
enabled by a small difference in the phase of the valley-orbit coupling between the two dots, and it 
in turn allows controllable rotations of two-dot valley eigenstates at a level anticrossing. We present 
a comprehensive analytical discussion of this process, with estimates for realistic structures. 



I. INTRODUCTION 

Since its inception three decades ago, quantum com- 
puting (QC) has developed energetically and spurred in- 
cessant innovation, uniting researchers from various areas 
of physics One area that has witnessed steady progress 
in recent years has been quantum control of systems 
at the nanometer scale, where state of the art exper- 
iments have made manipulation of two-level quantum- 
mechanical systems a reality. At the same time, in 
addition to accurate and reliable control, a quantum 
computer requires scalability and long coherence times. 
This requirement has highlighted solid-state spin sys- 
tems as a natural choice for quantum bits (qubits)pHSl 
Within this area, Si has been known for its extraordi- 
narily long coherence times thanks to the absence 
of piezoelectri c ele ctron-phonon coupling, 1 - 1 weak spin- 
orbit coupling 14 ! 15 ^, and nuclear-spin free isotopes, al- 
lowing removal of the hyperfine interaction by isotopic 
purification!^ As a result, Si spin QC has emerged as 
an active subfield of modern condensed matter physics. 
Outstanding experimental progress in Si spin QC has 
been reported in the last few years in Si quantum dots 
(QDs)pH2H and in donor-based architecturesPEU The- 
oretical research on Si QDs has also evolved at a brisk 
p ace IHH6i] Concomitantly, QDs in other group IV ele- 
ments are being explored for QC: carbonp2H£3 including 
nitrogen- vacancy centers in diamondpSl and GeP^ 

Group IV materials (C, Si, Ge) are notable for having 
equivalent conduction band minima known as valleys. In 
low-dimensional Si nanostructures typically the two low- 
lying valleys perpendicular to the interface are impor- 
tant. The valley degree of freedom, investigated many 
years agop^ has received renewed attention in the con- 
text of QC. The multiplicity of the Hilbcrt space brought 
about by the existence of equiv alent valleys has been 
shown to hamper spin QCPHEH At the same time, the 
interface potential gives rise to a valley-orbit coupling, 
which has been studied extensively in recent years, both 
experimentalljEsHUl and theoretically!^!!!] 

Addressing specific valley eigenstates is a profound, 
challenging and unresolved problem. Despite appar- 
ent similarities, valley physics is significantly different 
from spin- 1/2 physics, and a pseudospin picture of val- 
ley physics is of limited utility. Firstly, no overall con- 



servation law exists for valley composition as for spin. 
Secondly, a comparison of valley-orbit coupling and spin- 
orbit coupling sheds additional light on this distinction. 
In the absence of spin-orbit interactions the one-electron 
wave-function factorizes into an orbital part and a spin 
part, whereas no such separation exists for the valley de- 
gree of freedom. Unlike the spin-orbit interaction, the 
valley-orbit coupling cannot be viewed as a wave vector- 
dependent effective Zeeman field that can tune valley 
dynamics. Finally, even if the valley splitting is large 
enough to be resolved, no clear unambiguous signature 
of individual valleys exists. A smoking-gun experiment 
to identify valley-split states is an ambitious target, and 
several checks must be made on sets of quantum mechan- 
ical states to prove that they are indeed valley-split pairs. 
For example, for one electron in a single QD it was shown 
that relaxation between different valley eigenstates pro- 
ceeds on longer time scales.^ For two electrons in a single 
QD the ground state moves down as a function of mag- 
netic field. EHI and in a double QD (DQD) different valley 
eigenstates may be identified via resonant tunneling.^ 

Since the two valleys are separated by a wave vector of 
the size of the Brillouin zone, manipulation of the valley 
degree of freedom is a difficult task. No scheme has been 
experimentally demonstrated to date for achieving rota- 
tions of valley eigenstates in a 2DEG or in a single QD. 
Nevertheless, a recent work demonstrated that rotations 
of valley eigenstates can be implemented in a DQDP^in 
which one can engineer local variations in the magnitude 
of the valley-orbit coupling. This is due to the fact that 
the magnitude of the valley-orbit coupling in Si can be 
enhanced by a gate electric field P*l 

The valley-orbit coupling is a complex number, with 
a magnitude and a phase. The meaning of the phase is 
somewhat elusive, and it is not observable in a single QD 
- one requires at least a DQD to observe phase-related 
effects. The gate electric field affects both the magnitude 
and the phase of the valley-orbit coupling, even though 
its effect on the phase is much smaller.^ In this work 
we demonstrate that controllable valley rotations in a 
DQD can also be accomplished using the phase of the 
valley-orbit coupling. We focus on a different parameter 
regime from that considered in Ref. 1921 One key con- 
cept that we exploit is the fact that tunneling between 
like valley eigenstates (i.e. conserving the valley eigen- 
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state index) and between opposite valley eigenstates (not 
conserving the valley-eigenstate index) are not indepen- 
dent processes. Rather, one occurs at the expense of the 
other, and by controlling the phase of the valley-orbit 
coupling one can tune the tunneling between these two 
regimes. Control of the phase of the valley-orbit coupling 
can be achieved using a gate electric field. We devise an 
analytical model for the valley-orbit coupling, including 
the correction due to the gate electric field. We deter- 
mine expressions for the tunneling parameters between 
valley eigenstates, and their dependence on the phases 
of valley-orbit couplings in the two dots, as well as nu- 
merical estimates and a feasibility study for experimental 
implementation. We find that the rotation most likely to 
succeed experimentally involves the lowest two polarized 
triplet states in a DQD, which have different valley eigen- 
state composition. 

The outline of this paper is as follows. In Section |TT] 
we review briefly the model of the DQD, presenting the 
confinement potential and envelope functions, as well as 
an analytical scheme for calculating the valley-orbit cou- 
pling. In Section |III| we focus on the tunneling between 
valley eigenstates, determining the effect of a gate elec- 
tric field on the phase of the valley-orbit coupling and 
its subsequent effect on interdot tunneling. We supply 
numerical estimates of the change in phase of the valley- 
orbit coupling and of the intervalley tunneling parame- 
ter. Section |IV| discusses the coherent rotation of valley 
eigenstates of the two lowest polarized triplets in a DQD, 
including experimental considerations. We conclude with 
a summary and outlook in Section [V] 



II. SINGLE AND DOUBLE QUANTUM DOTS 

For consistency, we provide in this section the de- 
tails of m ultivalle y quantum dots, following previous 
discussions! 75 * 91 * 9 -^ An electron in a single quantum dot 
D experiences the potential 



the EMA equation for mot ion in the z-direction, per- 
pendicular to the interface) 86 * 91 * is a variational wave 
function,^ 



ip(z) = M e~ 6{-z) + N (z + z ) e" 



(3) 



where fc/, = y 2m f ^ J " , kgi is a variational parameter, and 

continuity of ip[z) at z = requires M = Nzq. In the 
basis {D^}, the Hamiltonian Hp, describing one electron 
in one multivalley QD, has the form 



H 



D 



A D 
A* D 



(4) 



with the confinement energy Ed and valley-orbit coupling 
A D = |A D |e-^ D , and e D > \A D \. The valley-orbit 
coupling in dot D is given by^S 



Ad = (D z \{U 6(z) + eF D z)\D g ) 



The eigenstates of H D are 



\D+) = 



1 



(\D x )±e i 4">\D s )). 



(5) 



(6) 



This study will focus on a double quantum dot, with 
the left dot located at xl — ~xq and the right dot at 
x r = x o- Henceforth we use D = L,R quite generally 
unless we need to refer specifically to the L. R dots. We 
assume the presence of a top gate that can be adjusted 
independently for the L, R dots, thus Fjj has different 
values for the two dots, Fl and Fr respectively. Taking 
into account these features, the confinement potential for 
this DQD can be written in the form 



Vdqd = 7r~^ \ Min [0 ~ x o) 2 , {x + x ) 2 } + if 



(7) 



+ eEx + eF L z + eF R z + U 0(z), 



V D (x,y,z) 



2m* a 2 



(x - x D ) 2 + y 2 



-U a 9(z) + eF D z. 

(1) 

The location of the dot is given by (xd, 0, 0) and its Fock- 
Darwin radius is a, with m* the Si in-plane effective mass, 
Fd the interface electric field, and Uq 9{z) the interface 
potential, with the interface at z — and 9(z) the Hcav- 
iside function. In a multi-valley system, in the effective 
mass approximation (EMA) the QD wave functions are 



Dt{x,y,z) = (j> D (x,y)ip(z)u e {r)e 1 ' 



(2) 



where u^(r) is the lattice-periodic Bloch function, the 
valley index £ = {z, z}, and k z S — ±fc , with k = 
0.85 {2-K/asi) and asi the Si lattice constant. The z and 
z states have a vanishingly small overlap which is ne- 
glected. The envelopes 4>n(x,y) are Fock-Darwin states, 

_(x-x £ ,) 2 +h 2 

4>D{x,y) = [l/^av 7 ?)] e ^ . The solution ip{z) of 



where the in-plane electric field E gives the interdot de- 
tuning. The overlap of the single-electron wave functions 
/ = (L(\R() 7^ 0, which motivates us to construct or- 
thogonal single-electron wave-functions as in Refs. 151751 
These wave functions are denoted by \Lc) — Ji^Lj^kL 



and \R%) 



l ^ ) ~ g|Li> , where g = (1 - \/T^P)/l, so 
that (Lf|-Rf) = 0. Next, we orthogonalize \L±) and 



\R±) as \L ± ) 



\L±)-g\R±) 



and \R±) 



\R±)-g\L±) 



^/l-2lg+g2 ' ^/l-2lg+gi ' 

Henceforth we use the states \D±) and all matrix ele- 
ments carry a tilde (though matrix elements with and 
without tildes are numerically almost identical.) In this 
basis we have a slightly modified valley-orbit coupling 
A = (D z \Uq6(z) + eFz\L>z) as well as modified confine- 
ment energies Ed- The assumptions underlying this for- 
mulation of the problem have been discussed at length in 
Ref.HU 
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III. INTERDOT TUNNELING BETWEEN 
VALLEY EIGENSTATES 



neling between different valiey eigenstates which can be 
effective in the neighborhood of a levei anticrossing. 



Electrostatic interactions alter the valley-orbit cou- 
pling, as was shown in Ref. |HS1 Even though the frac- 
tional change in the valley-orbit coupling can be sub- 
stantial, electrostatic interactions cannot rotate different 
valley eigenstates into each other on a single QD. Mag- 
netic interactions are much slower than electrostatic ones 
due to the smallness of the Bohr magneton. Consider a 
linear magnetic field gradient [Ib{9B / dz) z. The form of 
the valley-orbit coupling matrix element due to this in- 
teraction is the same as the electrostatic potential eFz, 
except it is smaller by ~ 10 7 even for the l arges t magnetic 
field gradients achievable experimentally. 1 -^ Therefore, 
unsurprisingly, no scheme has been devised for manip- 
ulating valley eigenstates via control of the valley-orbit 
coupling in a single Si QD (but see Ref. [65] for C) . 

This work is primarily devoted to a new method of 
harnessing the (small) electric field effect on the phase of 
the valley-orbit coupling in order to achieve controllable 
rotations of valley eigenstates in a double quantum dot. 
When the top gate electric field is the same for the L 
and R dots, so that Fl = Fr, the effective confinement 
potential experienced by the two dots is the same. If, in 
addition, the interface is sharp along the growth direc- 
tion and flat perpendicular to it, or if interface roughness 
is correlated over distances much shorter than the size of 
the QD, the valley eigenstates \D±) are identical in both 
dots!^ Under such circumstances interdot tunneling oc- 
curs only between the same valley eigenstates (+ to + 
and — to — ), while interdot tunneling between valley 
eigenstates (+ to — and — to +) is suppressed. In this 
section we discuss the modification of the phase of the 
valley-orbit coupling due to a top-gate, and demonstrate 
that tuning the top-gate electric field to be different on 
the L and R dots enables a small amount of interdot tun- 



A. Gate effect on valley-orbit coupling 

To determine the valley-orbit coupling A^ of Eq. 
we expand the lattice-periodic functions u^(r) as 



iK- 



(8) 



K 



with K reciprocal lattice vectors. For a perfectly smooth 
and perfectly sharp interface, neglecting the effect of Fd 
for the time being, A^ will be referred to as the global 
valley-orbit coupling Aq, which can be expressed as 



A = U N 2 z 2 



K,Q Z 



£* -5 
"k c k+q, 

h + iq z 



(9) 



where q z = Q z — 2ko. For practical purposes, since the 
Umklapp terms do not contribute, Ao can be accurately 
approximated by the term with Q z — only 



An 



U,N 2 z 2 
fct, — 2iko 



K 



K< K 



(10) 



We note that Zq depends on kgi, which in turn depends 
on the gate electric field. Nevertheless, the effect of the 
gate on zq is weak, and zq is mostly determined by kf,, 
which is fixed by the interface potential. 

Our first task is to determine the effect of the top gate 
Fd on the valley-orbit coupling. We evaluate the contri- 
bution due to Fd, which arises from the matrix element 



A B = {D z \eFz\D_ z 



/>OG nOO 

<F ! I dxdy\(f)(x,y)\ 2 I dz z\^(z)\ 2 e~ 2lk ^ z u* z {r) u_ z {r) 



— oo J — oo 



(11) 



eFN 2 



z 2 dz z eV° b+iq ' )z + / dzz{z + z ) 2 e-^- iq ^ z \. 



The integrals are trivial (see Appendix IB]) , and terms with Q z ^ are negligible, yielding 



Az^eFN 2 ^^ 

K 



+ 



4z 



(k b + iq z ) 2 (k Si -iq z ) 2 (k Si ~ iq z ) 3 (k Si - iq z ) 4 



The full valley-orbit coupling, including the gate correction, can be written as 



A D « N 2 (Y.CkCk){ 



Uoz, 



oz 



kb — 2iko 



eF 



4zn 



{h + iqz) 2 (ksi-iqz) 2 (ksi-iq z ) 3 (ksi-iqz) 4 



(12) 



(13) 
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The phase <f>r> 



arg Ad of the valley-orbit coupling 
Im Ad 



On = arctan 



Re A, 



(14) 



Assuming that the only difference in between the 
dots comes from the gate electric held over one dot (which 
we take to be L), so that = Aq + A e while A# = Ao, 
one can write approximately 



i>n = arctan 



arctan 



/ImAi 
VReA L 

Im (A 



arctan 



/ImAj 
\Re Ai 



A, 



Re(A + A E ) 



arctan 



ImA 



,ReA 
(15) 

Experimentally, both top and back gates are required to 
control the valley-orbit coupling and QD energy levels 
independently, as discussed extensively in Ref. l92l 




Figure 1: Interdot tunneling between valley eigenstates for 
Si/SiC>2. The parameters used are mj = 0.4mo,(7o = 
3.0eV, a = 8.2nm and b = 1.060nm. 



B. Tunneling between like and opposite valley 
eigenstates in a DQD 

Two tunneling parameters are relevant to a multivalley 
DQD. Tunneling between like valley eigenstates is given 
by t , which we expect to be dominant g iven the out- 
standing quality of present-day Si interfaces! 20 * 26 ! Tunnel- 

ing between opposite valley eigenstates is given by t |_. 

This latter parameter is zero if the valley-orbit coupling 
is exactly the same on both dots. The general formulas 
for t and t l are 



* r 

= 2^ 



^-i(4>l-4>r) 



— - [1 — e -»(*i-^H)] 



(16) 



where t = to + s, with to = (L^\H\R^) and the two- 
particle term s = (L^ Lf^VeelL^ Rf 1 ), with the su- 
perscript (i) denoting the i-th electron. The tunneling 
matrix element t, including the Coulomb contribution is 
calculated in the Appendix. This tunneling parameter 
is the same for both values of £. The gate electric field 
Fd affects both the amplitude and phase of An, and by 
tuning Ajj to be different on the two dots the tunneling 

parameters t and t |_ can also be tuned using a gate. 

One can keep Fr constant and tune Fl so that <f>L is 
different from <f>n. Assuming this difference to be small, 
we expand the exponential e^ % ^ L ^^ R ' . The parameter 

t acquires a correction linear in the electric field, but 

it remains effectively t, and we shall assume t ss t 

henceforth. However, t |_ can be expressed as 



t- 



it ,~ 



arctan 



Im (A 



Re (A + A E ) 



arctan 



Im A 

ReA 
(17) 



Table I: Numerical estimates of Ao, Ae and t |_ in Si QDs. 



Material 


Si/SiGe 


Si/Si0 2 


m6=0.4mo 


mt=0.3mo 


\Ao\(peV) 


200 


230 


150 


\A E \(peV) 


1.7 x 10 -2 


6.9 x 10~ 2 


7.9 x 10~ 2 


\4>L — (f>R,\(rad) 


8.6 x 10" 5 


2.9 x 10" 5 


1.2 x 10" 4 


\i-+\(neV)(d = 2.5) 


0.28 


0.20 


0.84 


\t-+\(fj,s) 


2.3 


3.3 


0.8 



By varying the gate electric field on one of the two QDs, 
t |_ can be tuned from zero to the desired time scale. Ex- 
perimentally, t and t |_ can be identified using resonant 

tunneling as described in detail in Ref. 91. 

We discuss next a series of numerical estimates for the 
valley-orbit coupling, the effect of an electric field on it, 

and the interdot tunneling parameter t |_. The dimen- 

sionless quantity d — Xq/cl represents half the interdot 
separation in units of the QD radius. A realistic inter- 
facial electric field of 150 kV/cm is assumed to obtain 
numerical estimates. For a Si/SiC>2 interface, with Uq ~ 
3eV and mj, = 0.4too, b — 1/ksi is optimized at 1.06 nm, 
and we find |Ao| ~ 230/ieV. For a Si/SiGe interface, with 
Uq ~ 150meV and mj = m , b is optimized at 1.217rtm, 
and we find |A | 200/ieV". Table ffl displays A , A E and 
t |_ for experimentally realistic parameters, and the pa- 
rameter t |_ is plotted in Figure 111 The key effect investi- 
gated in this work concerns the phase change due to A^;, 
which enables tunneling between opposite valley eigen- 
states. Determination of t (_ for experimentally relevant 

parameters yields a time scale of 0.8 — 3.5/^s, which can 
be easily accessed in the laboratory. This demonstrates 
the feasibility of two-dot manipulation. 

We note in closing this section that the phase of the 
valley-orbit coupling cannot be of use in the case of 
a single QD, regardless of its occupation number. In 
fact the phase of the valley-orbit coupling cannot be 
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{&, fff , f f^, fff , Tf «} of Ref. [7S] 
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Figure 2: Triplet levels in a double quantum dot in the ab- 
sence of interface roughness, plotted as a function of the de- 
tuning, which has been suitably non-dimensionalized as in 
Ref. 1751 We have assumed a bare valley splitting Ao of 
O.fmeV. 



measured in a single QD. Both the confinement energy 
£o = (D Z \T + Vd\Dz) and the on-site Coulomb energy 
it are independent of the phase of A. To extract infor- 
mation about the phase of the valley-orbit coupling the 
minimal requirement is a DQD, in which the phase dif- 
ference between the two dots gives rise to the measurable 
quantity t |_, as outlined above. 
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Barring a constant offset and noting that t + _ 
this basis the Hamiltonian is 



/— 2|A| 

2|A| L+ 
H T = t 
-t 
t*_ + t-i -5 ) 



(18) 



t-+, in 



(19) 



The detuning is defined as is defined as 5 = (ex — er) — 
(u — k), where u is the on-site Coulomb repulsion and k 
the two-electron two-dot direct Coulomb integral, both 
discussed in Ref. [75] We diagonalize this Hamiltonian 

treating as a perturbation. The cigenstatcs without 

< + _ are f±± and 



IV. COHERENT ROTATIONS OF 
TWO-ELECTRON STATES IN A DOUBLE 
QUANTUM DOT 

Two-electron states may be spin singlets or spin 
triplets. Previous work has found that the two lowest 
energy spin-singlet levels do not cross as a function of 
detuningPS Even for the case in which the lowest sin- 
glet state is easily initialized, no value of the detuning 
exists where one can controllably induce mixing between 
this state and another valley eigenstate of different val- 
ley composition. Therefore we focus on the spin triplet 
branch, in which it is easiest to access the polarized spin- 
down triplet states by applying a magnetic field of or- 
der 1-2 T. When the lowest-energy spin-polarized triplet 
state is initialized, the detuning can be swept to an anti- 
crossing in the triplet energy-level spectrum, where one 
can mix two spin-polarized triplet states with different 
valley composition. This provides a scheme for the im- 
plementation of coherent rotations of valley eigenstates, 
as we will see below. 

With the magnetic field set to zero for the time be- 
ing, we investigate the dynamics of the lowest two triplet 
states. The difference in |A_d|, i.e. |A/,| — |A^|, induced 
by the gate electric field can be ignored for our purposes. 
We therefore write |Al| = |A^| = |A|. We use the basis 



T> = 



Tl 



e> 2 + 2t 2 V £ o 



t ~ 



T LR_ 



1 f 



LR 



T 



RR 



~<2 



t ~ 



T LR_ 



2t 2 ^ £ o + 



*-f 



LR 



T 



RR 



1 



rfisym _ ■*■ / rfiLR, , rfiLR 

v'2 ! ' 



The notation used above is 

-6 + VS 2 + 8t 2 



-S-VS 2 + 8t 2 



(20) 



(21) 



Reliable preparation of the state T_j5_ is accomplished 
by raising the detuning S so that the right dot is much 
lower in energy than the left. A magnetic field of 1-2 
T is applied to separate the spin-polarized triplets from 
the other two-electron states. At this point two electrons 
can be unambiguously initialized on the right dot. This 
has already been demonstrated experimentally.^ (When 
a magnetic field is applied, the orbital states do change, 
but the effect will be the same for all triplet states - 
there will be a constant offset, which does not affect the 
argument presented here.) 
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The energies of the spin triplet states as a function 
of detuning are plotted in Fig. [2j The lowest energy 
states are and T^. Referring to Fig. [5J one now 

sweeps the detuning until the two states and Tt- 
cross. At the point when they cross, T+_ ~ T+^. Thus 
the matrix element mixing the lowest two triplet states is 

effectively , which gives an anticrossing at this point. 

Since usually A S> t, the detuning at the anticrossing 
is approximately given by 5 « 2A (corresponding to the 
energy difference of the two triplet states when S <C 0). 
The tunneling matrix element between T_j5_ and at 
the anticrossing is 

(fL R \H T \f<_)= £ o (fj- fi |g T |fff) 
V^o <2 + 2? 2 

/ 5 + VP + 8t 2 \ - 

^ \jzP + 2~8\/P + 8^ + 16i 2 ' 

Since (5 ~ 2A 3> f, the above tunneling matrix clement 

is approximately t |_. The width of the anticrossing is 

2y2t h- At the anticrossing point the two states, T5 

and Tr^, can be rotated into each other on a time scale 

given by 7i/f [_, which constitutes a coherent rotation of 

the valley state of two electrons. The values of t |_ given 

in Table |TJ which are of the order of /its, can be regarded 
as maximum values for the given parameters. Tuning 

the gate electric fields Fl and Fr allows one to tune t |_ 

from zero to the maximum achievable value, controlling 
the time scale of the rotation. 

Experimentally one needs to know the valley-orbit cou- 
pling A in order to determine where the anticrossing oc- 
curs. For this scheme to work it is necessary for the 
valley-orbit coupling to be measured beforehand! 30 * 79 * In- 
terface roughness gives an additional electric field, but 
that is a static offset. Finally, applying a large gate elec- 
tric field may modify the number of electrons in the dot, 
thus this proposal requires a depletion dot, rather than 
an accumulation dot. 

The discussion above has focused on the implementa- 
tion of a coherent rotation of valley eigenstates, which in 
the language of quantum computing would constitute a 
a x gate. It is evident that the detuning <5 could provide a 
a z gate. It would appear that the method proposed here 
for rotating valley eigenstates has the advantage of longer 
pulses, which are better experimentally being more pre- 
cise. Nevertheless, appealing as it is, the scheme pre- 
sented in this work cannot be promoted to a qubit as it 
stands, the main concern in quantum computing applica- 
tions being dephasing. Longer pulses require a longer co- 
herence time T£ (as found in the context of the hyperfine 
interaction in GaAs.^J) We have assumed the hyperfine 
interaction to eliminated through isotopic purification, 
phonons to be frozen out at dilution refrigerator temper- 
atures, while interface roughness and screened charged 
impurities give rise to a constant offset in the QD confine- 
ment potential. This leaves the biggest source of dephas- 



ing as charge noise due to dangling bonds in the vicinity 
of the interface, which cause fluctuations in the detuning, 
and give rise to random telegraph and 1/f noise Near 
the anticrossing, where the coherent valley rotation is 
performed, the energy splitting between the two levels is 
approximately linear in 5, thus the system is susceptible 
to charge noise [this term oc S arises from matrix ele- 
ments of the form (L-\Vdqd\L-) and (R±\Vdqd\R±) , 
i.e. one of the states has occupation (1,1) whereas the 
other has occupation (0, 2).] The absence of a sweet spot 
in the qubit energy splitting does not allow one to mini- 
mize the sensitivity to noise. Therefore, unless this noise 
can be filtered out the coherence time TJ may be too 
short for quantum computation. Reliable numbers for 
1/f noise in single and double Si QDs are scarce, though 
1/f noise and the charge offset drift has been measured in 
SETsPand are indicative of the results to be expected in 
QDs. At the same time, there has been progress of late in 
combating the effect of noise in QD spin qubits, such as 
a singlet-triplet qubit via dynamical decoupling,^ com- 
posite pulsesp^ and by growing a buried quantum dotP^ 
further from the interface. 

A recent publication^ has demonstrated that quan- 
tum bits can be implemented using the valley degree of 
freedom in Si QDs. One common feature of this work and 
Ref. I5H is their reliance on local control of the top gate 
over one quantum dot in order to generate a difference in 
the valley-orbit coupling A^ for the two quantum dots in 
order to mix two polarized triplets with different valley 
composition. Nevertheless, the two schemes employ dif- 
ferent principles and are operated in different parameter 
regimes. In Ref. EH it is the difference in the magnitude 
| A# | of the valley-orbit coupling that is exploited in or- 
der to mix the two triplet states Tf 3 ^™ and T"™* 1 in the 
far detuned regime 5 <^ 0. In contrast, this work exploits 
the difference in the phase of Ac in order to mix T_^_? 
and Tt- at a point where they cross, and it is not A^ 
itself that causes the mixing but the tunneling between 
valley eigenstates that it gives rise to. The conclusion 
that emerges from these two works is the following: the 
electric field changes both the magnitude and the phase 
of A. When one uses the change in magnitude of A, the 
change in phase is irrelevant. When one uses the change 
in phase of A, the change in magnitude is irrelevant. 



V. SUMMARY AND OUTLOOK 

We have demonstrated that controllable rotations be- 
tween valley eigenstates can be performed in double 
quantum dots in Si using gate electric fields. Finite in- 
terdot tunneling between valley eigenstates is enabled by 
a small difference in the phase of the valley-orbit cou- 
pling in the two dots, and it in turn allows controllable 
rotations of two-dot valley eigenstates states. 

The scheme we propose has relied on a DQD. The pos- 
sibility exists that extensions may be devised for systems 
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of three QDs, which have also been studied recently!^^ 
Such extensions may uncover new and interesting physics 
and will be addressed in the future. Furthermore, the 
scheme we propose is expected to ha ve analogs in car- 
bon nanotube and graphene QD d 62 * 66 ^ as well as recently 
realized Ge QDsP 

This work is supported by LPS-NSA-CMTC and by 
the National Natural Science Foundation of China under 



grant number 91021019. We thank N. M. Zimmerman for 
pointing out the role of charge noise in this system. We 
gratefully acknowledge discussions with S. Das Sarma, 
Andre Saraiva, H. W. Jiang, Andras Palyi, Wang Yao, 
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Appendix A: Interdot tunneling parameter t 

The interdot tunneling parameter t = to + s. The single-particle part of the tunneling parameter, to, does not 
depend on the form of the z-wave function, and is the same as that calculated in our previous papers (Ref. I75p . We 
therefore require the Coulomb integrals for the tunneling parameter. We wish to evaluate the integral 

, n ™ 1 f rt f e i> } e ?» W 2l V(z 2 ) 2 tk , 

U(P,Q) = — / d 3 n / d 3 r 2 V{ V ^ 21 , Al 

it 2 a 4 J J \r 1 -r 2 \ 

where ip(z) is real. Define 

d 3 k 



F(k) = J d 3 re lkr f{r). 



The Fourier transform of /(r) = 1/r is F(k) = Aw/k 2 . We write the term 1/|t*i — r 2 \ as a Fourier expansion 

1 47T f q e -ik2-{ri-r 2 ) -y f ^ e -ik 2 -(r 1 -r 2 ) 

\, 2 ^ M J "<2 



We Fourier transform the x and y-dependent terms 



e * = - a a K 1X 



.2J 



j2 lr 

4-7T 

\2 



The Fourier transform of ip{z) will be called ip{k z ) for now, and evaluated explicitly later. 

1 f°° 

V^i) 2 = 7T / d*i»e- <fcl ' a V(*u) 



2tt 



2tt 

The Fourier transform ip(k z ) is given by 



— oo 



/0 /"oo 
dz e ( fe "+^)^ + N 2 dz(z + z Q ) 2 e -V>si-ik.)* 
-oo JO 

N 2 z 2 N 2 z 2 2N 2 z 2N 2 



(k b + ik z ) (k Sl -ik z ) (k Sl -ik z ) 2 (k Si -ik z ) 3 ' 

The integral I(P, Q) becomes 



I i>( R )= Mil p e* 1 ^* 31 ^ * e ^ ip(k lz )ip(k 3z ) 

(A7) 



(A2) 



In - r 2 | - (2.)3 y d3fc2 ^1 = 2^ 7 d3fc2 ■ (A3) 



(A4) 



(A5) 



(A6) 



1 /" rf3fc l r ik^-(l> Qi, 

2tt 2 / A; 2 



,2 t 2 



1 1 
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where R = P — Q. Since ip(z) is real, ip{ki z )ip(—ki z ) = iV^iz)! 2 - The Coulomb tunneling parameter in our setup is 



s = -imx x) k2 



e ik ± X cos#g- 



i t / m9 e 2 l f dk±dk z k± . ,, . ° 2fc l , ~ . . , 9 , , . 
IV>(M| 2 = — y fc 2 + fc2 Jo(fc±^o)e-^|^(fe,)| 2 , (A8) 



recalling that s « s. This integral can be evaluated numerically for specific values of fcs. 



Appendix B: As integrals 



The four integrals required for the evaluation of are 



,c 



dzze {kb+iq ' )z = zl 



zl [ dz ze~ (ks '- iq ' )z = 
Jo 



f 

Jo 



e {k h +iq z )z 

(k b + iq z ) 



e (k b +iq z )z 

/ dz- — — — 

-oo ( k b + IQz) 



[ksi - Hz) 2 



2z / dzz 2 e-^ ksi - iq ^ z = 2z \ 2 / dz 



-(ksi-iq z )z 



4z 



dzz s e~ (ksi ~ iq " )z = 



(k S t - iq z ) J (k Si - iq z ) 3 

6 



P OO 

. dzz 2 e~^'^ z = . . 

{k S i - nz) Jo ( k si - n z r 

i 



£ (k b +iq z )z 

(h + iq z ) 2 



(h + iq z ) 2 
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